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1.0  INTRODUCTION 


This  report  describes  progress  during  the  second  year  of  a  research 
program  on  the  Fundamental  Processes  In  Partially  Ionized  Plasmas  con¬ 
ducted  in  the  High  Temperature  Gasdynamics  Laboratory  at  Stanford  Uni¬ 
versity.  This  research  is  supported  by  a  grant  from  the  Air  Force  Office 
of  Scientific  Research  (AFOSR-83-0108)  and  is  conducted  under  the  direc¬ 
tion  of  Professors  Charles  H.  Kruger,  Morton  Mitchner,  and  Sidney  Self. 
Three  Ph.D.  candidates  are  currently  conducting  their  doctoral  research 
under  this  program. 

Several  space  power  and  propulsion  systems  of  potential  long-range 
interest  to  the  Air  Force  Involve  partially  ionized  plasmas.  Such 
systems  include  MPD  thrusters,  both  open  and  closed  cycle  MHD  power,  and 
thermionic  energy  conversion.  Although  the  specific  configurations,  the 
exact  operating  conditions,  and  which  of  the  competing  systems  will 
prove  to  be  most  useful  in  the  long  term  remain  to  be  established,  it  is 
important  at  this  time  to  provide  a  broad  fundamental  research  base  in 
support  of  development  activity.  In  particular,  there  are  a  number  of 
key  issues  regarding  the  properties  and  discharge  behavior  of  partially 
ionized  plasmas  and  the  interaction  of  discharges  with  fluid  dynamics 
that  need  to  be  understood  before  the  potential  and  limitations  of 
competing  systems  can  be  fully  evaluated.  In  addition,  it  is  important 
that  outstanding  young  applied  scientists  be  educated  in  these  areas. 

The  present  research  on  partially  ionized  plasmas  is  discipline 
rather  than  device  oriented  and  is  currently  focused  on  three  major 
areas : 

1.  Recombination  in  molecular  plasmas 

2.  Discharge  effects:  plasma  electrode  interaction 

3.  Interaction  of  discharges  and  fluid  dynamics 

In  addition,  each  of  these  areas  involves  the  development  and  applica¬ 
tion  of  modern  plasma  diagnostic  techniques. 
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These  areas  are  overlapping  and  mutually  supportive.  Thus,  under¬ 
standing  of  plasma  properties  is  important  to  the  study  of  discharges 
and  their  interaction  with  fluid  dynamics.  In  the  same  spirit,  we  are 
interested  in  the  development  of  plasma  diagnostics  so  that  they  can  be 
applied  in  our  research  objectives.  In  each  area,  the  research  is  pri¬ 
marily  experimental  in  nature,  with  supporting  theoretical  studies  for 
the  planning  of  the  research  and  interpretation  of  the  data. 

Progress  during  the  second  year  in  each  of  the  three  research  areas 
is  described  by  means  of  Project  Summaries  in  Section  2.0. 

Publications  and  Presentations  resulting  from  this  work  are  cited  in 
Section  4.0,  and  Section  5.0  lists  the  personnel  who  have  contributed  to 
this  report. 


-  2  - 


Included  in  this  section  are  summaries  o£  progress  in  each  of  three 
project  areas.  Each  project  summary  contains  the  following  subsections: 
(a)  Introduction;  (b)  Research  Objectives;  and  (c)  Status  of  Research. 
Additional  descriptions  may  be  found  in  the  publications  listed  in 
Section  4.0. 

2.1  Recombination  in  Molecular  Plasmas 
Introduction 

The  goal  of  this  study  is  to  make  measurements  of  the  3-body  elec¬ 
tron  recombination  coefficient  a^T)  lor  the  process  e  +  M++  X^-*-  M  +  X^, 
and  to  compare  experimental  data  with  recent  theoretical  results.  Here  e 
denotes  an  electron,  M  and  M+  denote  an  atom  and  an  ion  respectively  of 
an  alkali  metal,  and  X^  is  a  molecular  species.  In  particular,  initial 
considerations  have  been  focused  on  the  alkali  metal  Cs  and  on  the  mole¬ 
cule  N2«  This  work  is  motivated  in  part  by  measurements  of  appreciable 
electron-density  nonequilibrium  in  plasma  boundary  layers  that  cannot  be 
adequately  explained  by  existing  data  for  known  recombination  processes 
Ll.l],  Such  boundary  layers  can  exist  in  devices  that  have  been  proposed 
for  generating  electrical  power  in  space  (and  could  exist  also  in  some 
space  propulsion  applications). 

The  recombination  of  an  electron  e  with  an  alkali  metal  ion  M+  by 
one  or  both  of  the  two  3- body  processes 

e+M++e*M+e  (1.1) 

e+M++Xk+M+e  (1.2) 

plays  an  important  role  in  many  applications  that  involve  the  nonequi¬ 
librium  behavior  of  partially  ionized  gases.  In  the  first  process,  the 
third  body  (which  is  required  in  order  to  conserve  momentum  and  energy) 
is  an  electron;  in  the  second  process,  the  third  body  is  an  atom  or  a 
molecule  X^  .  In  this  work,  the  temperature  dependent  recombination  rate 
coefficients  for  the  two  processes,  ae(T)  and  a^(T),  are  defined  by  the 
relations 


dt 


"VkVi 


(1.4) 
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where  n^  is  the  number  density  (cm  J)  of  ions  and  t  is  time  (sec). 

Experimental  studies  of  the  recombination  process  in  which  the  third 
body  is  an  atom  or  molecule  have  been  carried  out  using  (hydrocarbon) 
flames  [1.2j,  where  the  products  of  combustion  contained  several  molecu¬ 
lar  species.  To  extract  information  about  the  individual  molecules  X j  , 
the  assumption  was  made  that  the  contributions  from  the  different  mole¬ 
cules  in  the  mixture  were  additive.  Expressed  equivalently,  in  terms  of 
recombination  coefficients,  it  was  assumed  that  the  overall  rate  a 
(which  was  measured,  in  effect)  could  be  written  in  terms  of  the  indi¬ 
vidual  rates  Oj  by  the  formula 


*  *  1 XJ  °J  • 

where  Xj  ■  Oj/n  denotes  the  mole  fraction  of  species  Xj.  Here  nj  is  the 
number  density  of  species  Xj  and  n  is  the  total  number  density. 

The  results  obtained  in  this  fashion  have  been  strongly  criticized 
by  Bates  [1.3],  who  claimed  that  the  recombination  coefficient  for  a 
molecular  mixture  is  not  additive,  but  that  it  is  given  by  the  formula 
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The  quantities  Bj(p)  are  related  to  the  rates  of  de-excitation  of  M  by 
Xj  from  all  energy  levels  above  and  including  p,  to  all  levels  below  p. 


The  factor  n  is  given  in  terras  of  atomic  constants  and  the  temperature  T 


as  (2TnnekT/h2)3/2  -  KT3/2. 


For  a  single  species,  one  obtains  from  this  expresison  for  a  the 


result 


It  is  therefore  apparent  that  the  value  of  a  for  a  mixture  of  molecules 
cannot  be  written,  in  general,  in  terms  of  the  individual  a j •  To  be  able 
to  calculate  a  for  a  mixture,  the  factors  Bj(p)  must  be  used.  However 
before  being  able  to  rely  on  the  calculated  factors  Bj(p),  it  is  neces¬ 
sary  to  test  the  theory  by  obtaining  data  in  which  only  a  single  molecu¬ 
lar  species  Xj  is  present  in  the  recombining  gas  together  with  the 
alkali  metal  M.  It  is  clear  that  flame  experiments  are  unable  to  yield 
such  data. 


Besearch  Objectives 


1.  To  design  and  construct  an  experimental  system  that  can  produce  a 
recombining  alkali  metal  plasma  under  conditions  where  the  molecu¬ 
lar  three-body  process  dominates. 

2.  To  examine,  in  particular,  the  recombination  of  cesium  ions  in  the 
presence  of  nitrogen. 

3.  To  measure  the  desired  recombination  rate  as  a  function  of  the 
background  gas  temperature. 

4.  To  measure  the  recombination  rate  as  a  function  of  various  values 
of  the  partial  pressures  of  Cs  and  N£« 

5.  To  examine  the  recombination  rate  of  interest  in  a  mixture  of 
molecular  species. 

6.  To  design  and  construct  a  suitable  pulse  forming  circuit  for  a  high 
intensity  flashlamp  and  to  explore  the  use  of  pulsed  photoioniza¬ 
tion  as  a  means  for  producing  the  desired  plasma  conditions. 

7.  To  determine  the  electron  number  density  as  a  function  of  time  in  a 
plasma  afterglow,  and  to  deduce  the  recombination  rate  from  this 
measurement. 
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8.  To  develop  a  four  pin  probe  for  measuring  the  time  resolved  plasma 
electrical  conductivity  as  a  diagnostic  for  determining  the  elec¬ 
tron  concentration. 

9.  To  develop  a  method  employing  an  equilibrium  plasma  for  calibrating 
the  electron  conductivity  probe. 

10.  To  examine  the  possible  use  of  optical  diagnostic  methods  such  as 
measurement  of  the  emission  of  free-bound  continuum  radiation,  as  a 
means  for  obtaining  a  check  on  the  data  via  an  independent,  redun¬ 
dant  proceedure. 

11.  To  assess  the  possible  roles  of  the  present  of  the  dimer  ions  Cs2+ 
and  ion  clusters  such  as  N2Cs+  and  t^OCs*  'the  water  vapor  will  be 
present  in  small  quantities  as  an  impurity)  in  contributing  to  the 
measured  recombination  rates. 

Status  of  Research 

Research  during  the  period  February  1984  through  January  1985  pro¬ 
gressed  from  the  numerical  analysis  of  a  proposed  steady-state  super¬ 
sonic  nozzle  expansion  experiment  employing  a  high  temperature  plasma, 
to  conception  and  adoption  of  an  alternative  transient  photoionization 
afterglow  experiment.  During  this  period  the  analysis  of  the  afterglow 
experiment  underwent  extensive  development  and  much  of  the  necessary 
hardware  was  specified  and  procured.  This  status  report  presents  the 
rational  motivation,  fundamental  concepts,  and  major  hardware  designs 
which  were  produced  in  the  last  contract  year.  Included  herein  are  the 
analyses  by  which  the  limitations  of  the  nozzle  expansion  experiment  and 
the  advantages  of  the  pulsed  experiment  were  determined.  The  afterglow 
experiment  and  its  essential  hardware  are  described  in  detail. 

In  the  Annual  Report  of  March  1984  an  experiment  was  suggested  in 
which  an  arc  jet  would  be  used  to  produce  a  plasma  consisting  of  N2  (or 
some  other  molecular  species)  and  thermally  ionized  Cs  and  then  a  rapid 
expansion  using  a  supersonic  nozzle  would  be  employed  to  produce  steady- 


state  recombining  conditions.  Measurements  of  the  electron  number  den¬ 
sity,  ng,  before  and  after  the  expansion  would  be  compared  with  calcu¬ 
lated  values  obtained  from  the  electron  continuity  equation.  Since  the 
calculated  values  of  ng  depend  on  the  recombination  coefficient  (T) , 
such  an  experiment  would  provide  a  means  for  inferring  the  recombination 
rate. 


Two  non-dimensional  parameters  fundamental  to  any  recombination 
experiment  were  identified  in  the  last  Annual  Report.  One  parameter  was 
the  ratio  of  rates  of  recombination  on  molecules  to  recombination  on 
electrons,  r  =  aN2nN2^aene^  *  or<*er  tliat  ttie  recombination  be  domi¬ 
nated  by  the  molecular  third  body  process,  it  is  desirable  to  have  r  > 
10.  The  other  parameter  was  the  ratio  of  the  characteristic  time  for 
recombination  on  ^  to  the  residence  time  of  a  differential  volume 
element  of  plasma,  t  ■  (°tN2IlN2ne,:res^  ^ *  Since  the  equilibrium  and 
frozen  limits  correspond  to  t«1  and  t»1  respectively,  a  desirable 
value  for  t  would  be  t  =  1. 

Shown  in  Figures  1.1  and  1.2  are  the  results  of  numerical  calcula¬ 
tions  for  the  variations  of  electron  number  density  through  a  Mach  2 
converging-diverging  nozzle  for  various  representative  values  of  the 
upstream  stagnation  conditions  that  may  be  achievable  with  the  steady- 
state  gas  flow  system  presently  available  in  the  High  Temperature  Gas- 
dynamics  Laboratory.  The  operating  characteristics  of  the  arc  jet  heater 
in  that  system  limit  the  N2  partial  pressure  to  about  four  atmospheres, 
or  in  equivalent  terms  n^  <  10^  m-^  for  a  plenum  temperature  of 
3000  °K..  (Figures  lb  and  2b  are  for  conditions  which  are  somewhat  more 
extreme.)  The  term  effect”  refers  to  the  value  at  the  nozzle  exit 

of  the  fractional  decrease  in  ng  attributable  to  3-body  recombination 

with  N9,  that  is  In  (aM  =0)  -  n.j/n  .  To  obtain  these  results  the 
£  6  *12  6  6 

Hinnov-  Hirschberg  expression  [1.4]  was  used  for  ae(T),  and  the  theoret¬ 
ical  expression  proposed  by  Bates  et  al.  [1.5]  has  been  used  for  Ojj^(T). 
Sonic  flow  throughout  a  10  cm  long  nozzle  has  been  assumed  in  calcu¬ 


lating  residence  time  as  t. 


0.005//T. 


In  the  current  grant  period  the  necessity  of  determining  aN^(T) 
in  an  experiment  with  only  one  molecular  third  body  present  and  the 
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significance  of  the  parameters  r  and  t  has  remained  unchanged.  The 
gasdynamic  experiment  has,  however,  been  re-evaluated.  In  the  following 
sections  some  limitations  inherent  to  the  steady-state  expansion 
experiment  are  discussed.  As  an  alternative  approach,  a  pulsed 
experiment  which  creates  an  isothermal  plasma  via  photoionization  is 
shown  to  have  many  advantages. 

Limitations  of  the  Steady  State  Experiment 


It  appears  from  Figures  1.1  and  1.2  that  in  order  for  the  value  of 
ne  at  the  exit  of  the  nozzle  to  depart  significantly  from  the  frozen 
value,  it  is  necessary  to  have  t<1,  but  this  requirement  then  precludes 
achieving  the  condition  r>10.  If  ng  is  decreased  so  as  to  increase  r, 
the  value  of  t  is  also  increased,  thereby  reducing  the  effects  of  recom¬ 
bination  relative  to  frozen  flow.  This  behavior  may  be  seen  more  clearly 
from  the  relation 


(oN2aN2)  cres 


(1.5) 


obtained  by  combining  the  expressions  for  r  and  t  so  as  to  eliminate  ng. 
For  fixed  N2  partial  pressure  and  nozzle  dimensions,  the  right  hand  side 
of  (1.5)  is  a  function  only  of  the  stagnation  temperature.  Taking  PN  »4 
atmospheres  and  tres  ■  0.005//T,  one  obtains  r/x  ■  5.4  for  T  ■  2000  °K. 
and  r/x  *  7.2  for  T  *  3000  °K. 

Another  limitation  of  the  rapid  expansion  experiment  stems  from  the 
practical  difficulties  of  providing  optical  access  for  diagnostic  meas¬ 
urements  of  the  expanding  gas  at  more  than  a  few  axial  locations  (e.g. 
the  plenum,  the  exit,  and  perhaps  one  intermediate  position).  Since  the 
gas  temperature  decreases  by  approximately  1000  ®K.  from  the  upstream 
plenum  to  the  downstream  exit,  the  accuracy  of  determining  the  temper¬ 
ature  dependence  of  (T)  would  be  limited  by  the  relatively  coarse 
resolution  with  which  the  dependence  on  axial  location  of  ne  could  be 
measured. 


Rather  than  following  this  path,  an  entirely  different  approach  has 
been  examined  in  which  the  required  ionizational  nonequilibrium  is 
achieved  by  pulsed  photoionization  of  Cs  vapor  instead  of  by  a  rapid 
expansion  of  a  thermal  plasma.  This  alternative  approach  appears  to 
offer  significant  advantages. 

Description  of  Pulsed  Photoionization  Experiment 

The  radiation  source  chosen  for  photoionization  is  a  pulsed  short 
arc  Xe  lamp  directed  along  the  axis  of  a  tube  containing  the  test  gas 
(e.g.  Cs  vapor  and  N2).  A  schematic  of  the  arrangement  is  shown  in 
Figure  1.3.  The  Cs  vapor  and  N2  are  contained  in  a  one  inch  diameter, 
twelve  inch  long  inconel  tube  fitted  with  sapphire  windows.  The  Cs 
pressure  is  one  to  ten  torr  and  the  nitrogen  pressure  is  about  one 
atmosphere.  An  electric  tube  oven  heats  the  Cs-^  mixture  in  the  temper¬ 
ature  range  400  to  1200  X.  The  sapphire  window  has  75%  transmission  of 
radiation  down  to  200  nm  wavelengths  so  it  is  transparent  to  photons 
capable  of  singly  ionizing  ground  state  Cs  atoms,  i.e.  X  <  318  nm,  and 
opaque  to  photons  capable  of  ionizing  N2,  i.e.  X  <  200  nm.  The  Xe  flash- 
lamp  provides  a  2  ms  pulse  of  ionizing  photons  which  is  coupled  to  the 
Cs  test  cell  through  an  ellipsoidal  reflector. 

A  nonequilibrium,  isothermal  Cs  plasma  is  produced  by  the  radiation 
pulse.  The  free  electrons  are  rapidly  thermalized  through  collisions 
with  the  ^2*  As  the  electrons  begin  to  recombine  with  Cs*  in  the 
presence  of  the  nitrogen,  the  electron  number  density  would  be  deter¬ 
mined  as  a  function  of  time.  Possible  methods  for  obtaining  the  electron 
density  include  measurements  of  the  spectral  emission  from  free-bound 
and  high  lying  bound-bound  transitions  in  Cs,  as  well  as  measurements  of 
the  electrical  conductivity  using  a  four  pin  conductivity  probe.  The 
time  history  of  the  electron  number  density  is  directly  related  to  the 
rate  coefficient  for  recombination  on  N2. 


Experimental  conditions  would  be  chosen  so  that  the  recombination  on 
N2  proceeds  at  a  much  greater  rate  than  electron  loss  by  recombination 
on  electrons,  or  by  diffusion  of  electrons  to  the  test  cell  walls.  Under 


such  conditions  the  electron  kinetics  are  simply  described  by  the  dif 
ferential  equation 


■au  n  n  , 
N2  N2  e  * 


(1.6) 


which  has  the  solution 


(1.7) 


It  is  seen  that  the  slope  of  the  plot  of  l/nft  vs  time  gives  a«  (T) 
multiplied  by  the  known  concentration  of  N^.  The  temperature  dependence 
of  aN  (T)  is  empirically  determined  by  performing  the  experiment  over 
the  range  of  temperatures  that  can  be  achieved  with  the  apparatus  pre¬ 
viously  described. 

The  electron  number  densities  corresponding  to  the  parameter  value 
r  ■  10  are  given  in  Table  1.1.  Also  shown  in  the  table  are  character¬ 
istic  times  for  recombination  on  N2,  t^  ~  (a^  ny  ne)  * ,  for  electron 

2  2  2 
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diffusion  to  the  walls,  tdiff  ~  R  /°a  “  KR  n^/TT,  for  the  ionizing 
light  pulse  to  decay  from  a  peak  to  a  minimum  intensity,  tplg,  and  for 
thermalizatlon  of  free  electrons  tr^x  *  ^ve,N2* (Znig/m^)]  •  Here  R  » 
1  cm  is  the  radius  of  the  test  cell,  is  the  ambipolar  diffusion 
coefficient,  K  is  a  constant  that  depends  on  atomic  properties  of  Cs+ 
and  N2  [1.6],  ve  *-8  average  e,^  collison  frequency  [1.7],  and 

10.(2me/mjj^)  is  the  average  fraction  of  electron  kinetic  energy  lost  in 
all  collisions  with  N2  [1.8],  A  typical  value  for  the  flashlamp  pulse 
forming  network  is  used  for  tp^g. 

Table  1.1  shows  that  the  characteristic  times  are  ordered 
as  Trlx<'<'Tpls<'TN2^^T8iff  *  an< *  therefore  that  the  experimental  conditions 
result  in  an  isothermal  plasma  dominated  by  recombination  on  N2.  The 
degree  of  excitation  of  N2  by  encounters  with  photons  or  energetic 
electrons  is  negligible.  The  N2  molecules  do  not  significantly  absorb 
photons  with  wavelengths  greater  than  200  nm  11.9J.  The  energy  corre¬ 
sponding  to  an  electron  liberated  from  the  Cs  ground  state  by  a  200  nm 
photon  is  2.1  eV.  To  excite  a  ground  state  N2  to  its  first  vibrationally 


excited  state  requires  approximately  0.22  eV  [1.10],  or  each  electron 
can  excite  at  most  10  ^  molecules.  An  estimate  of  the  maximum  fraction 
of  N2  molecules  which  could  be  excited  in  this  manner  is  therefore 
10  •  <  ‘O’6. 


TABLE  1.1 

Characterization  of  Pulsed  Experiment 


Temp,  ‘K  ne>  cm"3  t^,  sec  tdlff,  sec  tpls>  sec  tflx>  sec 


400 

3 

X 

1010 

7 

X 

10"3 

0.9 

1 

X 

10“4 

4 

X 

10~9 

600 

7 

X 

1010 

8 

X 

10"3 

0.5 

1 

X 

10“4 

4 

X 

CT* 

1 

O 

800 

2 

X 

1011 

7 

X 

10"  3 

0.3 

1 

X 

10-4 

4 

X 

O 

1 

O 

Advantages  of  the  Pulsed  Experiment 


In  essence,  what  is  needed  to  determine  aM  (T)  is  a  measurement  of 
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the  recombination  rate  in  a  plasma  which  is  out  of  ionizational  equi¬ 
librium  and  which  contains  only  a  single  molecular  species  that  can  act 
as  a  third  body  for  recombination.  The  pulsed  photoionization  experiment 
satisfies  these  requirements  while  avoiding  many  of  the  technical  diffi¬ 
culties  inherent  in  the  steady-state,  rapid  expansion  experiment.  Since 
the  plasma  is  stationary  in  the  pulsed  experiment,  diffusion  (not  con¬ 
vection)  determines  the  residence  time  of  a  plasma  differential  volume 
element.  (Assuming  diffusion  is  the  major  cause  of  electron  loss  behind 
three  body  recombination  with  N2  implies  that  other  loss  mechanisms  can 
be  made  insignificant.  Loss  mechanisms  including  electron  attatchment 
and  cluster  formation  like  ^Cs*,  which  can  be  a  precursor  to  dissocia¬ 
tive  recombination  of  Cs2+»  are  currently  being  evaluated.  It  is 
believed  that  conditions  can  be  chosen  so  as  to  minimize  the  effects  of 
all  electron  loss  mechanisms  as  compared  with  diffusion  and  three  body 
recombination  on  N2.) 


The  characteristic  time  for  the  ambipolar  diffusion  of  a  Cs  plasma 
in  the  one  inch  diameter  cell  containing  N2  at  a  pressure  of  one  atmo- 


V  /  *  / 


*yuttze  is  about  500  milliseconds,  at  600  #K,  or  4000  times  longer  than 
the  convection-limited  residence  time  in  a  10  cm  long  sonic  nozzle  at 
2000  °K.  If  the  previously  derived  formula  for  r/x  is  used  with  tre8 
replaced  by  one  obtains  t/t  -  700  at  T  -  400  #K  and  r/x  -  250  at 

T  »  800  °K.  For  the  photoionization  experiment,  it  is  necessary  to  have 
t  1*(°h  ne)”*/tdiff  1  ^  1  (in  contrast  to  the  condition  x  *  1  for  the 
rapid  expansion  experiment).  It  is  possible  to  have  this  condition  on  t 
satisfied  and  at  the  same  time  have  r  >  10  in  the  pulsed  experiment. 

The  photoionization  experiment  creates  an  isothermal  plasma  because 

the  kinetic  energy  of  free  electrons  after  photoionization  is  quickly 

thermalized.  Recombination  occurs  at  a  constant  temperature,  and  the 

temperature  dependence  of  aM  (T)  is  measured  directly  by  simply 
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repeating  the  experiment  at  different  temperatures.  The  large  r/x  ratio 
obtainable  in  the  pulsed  experiment  allows  both  r  and  x  to  be  chosen  to 
simultaneously  eliminate  ambiguity  about  the  role  of  the  competing 
recombination  on  electrons  and  about  the  role  of  diffusion.  In  these 
respects  the  pulsed  photoionization  experiment  appears  to  be  more  favor¬ 
able  than  the  steady-state  rapid  expansion  experiment. 

The  Cesium  Cell 

Inconel  601  (Huntington  Alloys)  was  chosen  as  the  Cs  tube  wall 
material  because  of  its  exceptional  ability  to  resist  oxidation  and 
spalling  at  temperatures  up  to  1530  K.  The  1  inch  diameter  tube  will  be 
heated  at  atmospheric  pressure  in  a  tube  oven  to  temperatures  in  the 
range  between  400  and  1300  K,  and  it  must  withstand  oxidation  by  air 
externally,  and  reduction  by  Cs  Internally.  Where  the  inconel  tube 
extends  beyond  the  oven,  the  temperature  will  be  maintained  at  700  K  by 
two  Independent  ovens.  This  end  temperature  is  the  maximum  working 
temperature  of  the  sapphire  windows  (necessay  for  ionizing  and  diag¬ 
nostic  radiation),  instrumentation  feed-through  (for  the  conductivity 
probe),  and  vacuum  flanges.  Because  this  hardware  can  not  sustain  the 
entire  temperature  range  of  the  experiment,  the  cesium  cell  has  end 
regions  held  at  a  different  temperature  than  the  center  region. 


The  Cs  vapor  pressure  will  be  controlled  by  carefully  monitoring  the 
temperature  of  the  Cs  reservoir.  The  reservoir  will  be  maintained  20  to 
50  K  cooler  than  any  other  part  of  the  system  to  prevent  Cs  condensation 
in  other  parts  of  the  system.  The  oxygen-free  copper  reservoir  and  the 
N2  supply  and  vacuum  pump  will  be  attached  to  the  Cs  cell  by  heated 
stainless  steel  tubing.  The  diffusion  of  the  Cs  vapor  from  the  reservoir 
through  the  one  inch  diameter  by  12  inch  long  tube  to  90%  of  final 

concentration  will  occur  in  the  time  t»0.47  x  /D  [1.11].  Here  x  is  the 
tube  length  and  D  -  0.0911  (T/288)^2  cm2 /sec  is  the  N2-Cs  binary 

diffusion  coefficient.  For  T=*400  K,  t*1.85  hours.  Although  not  unaccept¬ 
ably  long,  this  mixing  time  will  be  decreased  by  natural  convection. 

The  N£  number  density  will  be  determined  from  simultaneous  temper¬ 
ature  and  pressure  measurements  before  the  cell  is  heated.  To  not  exceed 
the  1%  uncertainty  estimated  for  the  thermocouple  measurement,  a  1% 
limit  on  pressure  measurement  uncertainty  is  set.  This  accuracy  must 
hold  over  the  pressure  range  200  to  1500  torr  (resulting  from  a  final 

pressure  of  1500  torr  with  final  temperatures  ranging  from  400  to 

1300  K.)  These  requirements  will  be  met  by  a  series  of  Wallace  and 
Tieman  gages  which  span  the  desired  pressure  range. 

The  tube  and  reservoir  will  be  attached  via  a  bakeable  valve  to  a 
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liquid  nitrogen  trapped  diffusion  pump.  Pressures  of  10  torr  will  be 
achievable  with  this  system  in  order  to  remove  impurities  during  initial 
bake  out  and  subsequent  gas  changes. 

Photoionization 


Photoionization  of  the  Cs  by  Xe  flashlamp  radiation  occurs  through 
both  one  and  two  photon  processes.  The  single  photon  process  converts  a 
ground  state  6s  electron  to  a  free  electron  as  represented  by  the  reac¬ 
tion 


Cs(62S1/2) 


+  hv^  ♦ 


Cs+  + 


(1.8) 
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Here  the  photon  energy,  hvj,  must  exceed  the  ionization  potential  of  Cs, 
i.e.  A<318  nm.  Two-photon  ionization  first  excites  a  ground  state  elec¬ 
tron  with  a  resonant  photon;  a  second  photon  then  ionizes  the  excited 
atom.  The  two  steps  are 

cs(62S1/2)  +  hv2  C8*^2p3/2,i/2)  U.9) 

Cs  (62P3^2  Ay2)  +  hv3  ♦  Cs*  +  e  (1.10) 

The  resonance  wavelengths  are  12*852  nm  (J-3/2)  or  A2**894  nm  (J»l/2), 
and  ionization  from  the  6p  state  can  occur  with  any  photon  of  wavelength 
shorter  than  A3<508  nm  (J=3/2)  or  1^494  nm  (J*l/2).  The  cross  sections 
for  ionization  from  the  6s  and  6p  levels  are  approximately  10  cm  and 
10  cm  respectively.  Contributions  from  the  two  photon  process  to 
total  free  electron  production  are  estimated  to  be  substantial  [1.12]. 

The  experimental  arrangement  is  greatly  simplified  by  positioning 
the  light  source  away  from  the  Cs  and  the  heat  of  the  reaction  zone. 
Placing  the  source  outside  the  cell  necessitates  directing  the  light 
through  a  sapphire  window  at  the  end  of  the  Inconel  tube.  To  produce  a 
collimated,  intense  photon  flux,  a  high  radiance  lamp  is  required.  We 
have  chosen  a  five  atmosphere  Xenon  short  arc  lamp,  (L-2369,  manufac¬ 
tured  by  ILC  Technology  of  Sunnyvale,  CA.)  This  lamp  has  a  seven  mm  arc 
length  and  a  synthetic  fused  silica  envelope.  It  is  designed  for  pulsed 
operation  at  a  maximum  energy  input  of  300  Joules  per  pulse.  Pulsed 
operation  [1.13,  1.14]  and  the  synthetic  fused  silica  envelope  enhance 
ultraviolet  output  from  the  lamp.  The  short  arc  lamp  will  be  aligned 
along  the  axis  of  an  ellipsoidal,  aluminum  reflector.  (The  A-5000,  1000 
watt  reflector  is  sold  by  Photon  Technology  inc. ,  Princeton,  N.J.).  The 
arc  is  Imaged  at  the  far  focus  with  f/2.5  optics  by  positioning  the  arc 
at  the  near  focus  of  the  ellipsoid.  The  resulting  beam  diverges  with  a 
half  angle  of  11  degrees,  and  it  can  be  aimed  directly  through  the 
sapphire  window  into  the  tube.  Shallow  incidence  and  high  reflectivity 
of  the  polished  metal  tube  Interior  surface  will  insure  a  uniform  photon 
flux  through  the  cell's  length. 


A  multiple  mesh  pulse  forming  network  has  been  designed  to  provide  a 
roughly  rectangular  current  pulse  to  the  arc  lamp.  Pulse  width  will  be  3 
ms,  and  pulse  height  will  be  1,700  amps.  Begining  from  zero,  the  elec¬ 
tron  number  density  increases  toward  a  value  at  which  the  rate  of  crea¬ 
tion  of  electrons  by  photoionization  equals  the  loss  rate  due  to  recom¬ 
bination.  The  equilibriuum  electron  density  is  approached  asymptotically 
and  the  photon  flux  and  pulse  width  must  be  great  enough  to  produce 
sufficient  ionization  before  the  flash  terminates.  To  have  negligible 
interference  by  photoionization  during  the  measured  plasma  decay,  the 
current  pulse  will  fall  from  its  maximum  value  to  zero  in  0.9  ms. 
(Recombination  will  occur  on  the  time  scale  of  7  ms  ,see  tN^,  Table  1.1) 


Diagnostics 


A  four  pin  plasma  conductivity  probe,  developed  at  Stanford  in 
connection  with  research  on  MHD  plasmas,  appears  to  be  well  suited  for 
the  work  proposed  here.  The  principle  of  operation  of  the  probe  is  to 
measure  the  current  I  that  flows  when  an  ac  voltage  is  applied  between 
the  outer  pair  of  pins  (i.e.  electrodes)  of  a  linear  array  of  four  pins, 
and  to  measure  simultaneously  the  resulting  potential  difference  6V 
between  the  inner  pair  of  pins.  Since  J  (current  density)  is  propor¬ 
tional  to  1  and  6V  is  proportional  to  E  (electric  field) ,  it  follows 
from  Ohm's  Law,  E  *  J/o,  that  6V  ■  K  I/a  and  therefore  that  a  “  K 
dl/d(6V).  Here  a  denotes  the  electrical  conductivity  and  K  is  a  calibra¬ 
tion  constant  that  depends  only  on  the  geometric  features  of  the  probe. 
(In  practice,  6V  *  K.  I/a  +  constant,  but  this  constant  plays  no  role 
upon  differentiation.)  The  constant  K  can  be  determined  by  immersing  the 
probe  in  a  liquid  electrolyte  (e.g.  aqueous  KC1)  solution. 

The  conductivity  is  then  related  to  the  electron  number  density  by 
the  relation  a  m  (nge  )/(mevell),  where  veh  is  the  electron-heavy  colli¬ 
sion  frequency  and  e  is  the  electronic  charge.  Temporal  resolution  of 
50ys,  corresponding  to  operation  at  20  kHz  has  been  demonstrated 
[1.15].  Initially  this  probe  will  be  tested  on  a  steady  state  thermal 
plasma.  An  equilibrium  electron  concentration  of  3.6  x  10*®  cm~^  can  be 


achieved  by  heatiag  Che  Cs  cell  to  1300  K  and  the  Cs  reservoir  to  630  K 
corresponding  to  a  vapor  pressure  of  5  torr.  The  four  pin  probe  measure¬ 
ment  will  be  coapared  to  equilibrium  calculations  of  electron  concentra¬ 
tion  based  on  Cs  reservoir  and  cell  temperature  measurements.  The  probe 
will  also  be  coapared  to  optical  diagnostic  measurements. 

Measurements  of  the  recombination  radiation  continuum  have  been  used 
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to  measure  electron  densities  as  low  as  3X10  cm  in  the  afterglow  of 
a  Cs-Ar  discharge  [1.16,  1.17].  In  this  work  [1.16,  1.17]  a  monochro¬ 
mator  was  used  for  spectral  filtering.  By  employing  a  50  A  filter  in 
place  of  the  monochromator,  calculations  indicate  it  should  be  possible 
to  Increase  the  collected  radiation  by  a  factor  of  the  order  of  10  and 

therefore  to  lower  the  limit  of  measureable,  transient  ng  to  about  10^ 
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cm  .  This  limit  can  be  lowered  still  further  to  10  cm  when  measure¬ 
ments  are  made  on  a  steady  state  plasma.  Hence  optical  emission  diag¬ 
nostics  should  provide  a  check  for  the  conductivity  probe  during  testing 
on  the  thermal  plasma.  Upon  successfull  testing  in  a  steady  state 
plasma,  the  conductivity  probe  will  be  used  to  measure  electron  densi¬ 
ties  in  the  plasma  afterglow. 

At  present,  a  simulated  probe  and  high  impedance  amplifier  have  been 
built  and  tested  with  a  1  HI  simulated  plasma  resistance.  Successful 
performance  has  been  obtained  at  frequencies  up  to  500  Hz.  Higher  fre¬ 
quency  operation  is  hampered  by  capacitive  coupling  between  probe  leads. 
The  problem  of  stray  capacitance  is  currently  being  considered. 

Summary 


Because  of  the  limitations  of  the  nozzle  experiment,  effort  was 
directed  toward  analyzing  alternative  experimental  arrangements  for 
producing  the  required  nonequilibrium  conditions.  Among  the  alternatives 
considered,  the  approach  that  appeared  most  attractive  was  one  in  which 
the  lonizational  nonequilibrium  would  be  achieved  by  pulsed  photoioniza¬ 
tion  of  Cs  vapor  confined  with  ^  in  a  heated  test  cell.  The  recombina¬ 
tion  coefficient  aN  (T)  would  be  determined  by  measuring  the  decay  of  n 

tv  £  6 

as  a  function  of  time,  following  the  completion  of  the  photoionization 


pulse.  In  a  photoionization  experiment  the  recombination  is  non-steady 
In  time  and  uniform  in  space,  in  contrast  to  a  rapid  expansion  experi¬ 
ment  which  is  steady  in  time  and  nonuniform  in  space. 

Extensive  calculations  have  been  performed  that  demonstrate  the 
feasibility  of  the  photoionization  approach.  Possible  diagnostic  methods 
for  measuring  n£  in  the  range  of  plasma  conditions  expected  have  been 
reviewed,  and  a  method  based  on  measurement  of  the  electrical  conduc¬ 
tivity  (previously  used  in  this  laboratory  in  another  context)  appears 
feasible.  The  necessary  high  temprature  materials  plus  vacuum  and 
optical  components  have  been  designed  and  procured  and  6ome  diagnostic 
development  has  been  completed. 


Introduction 


The  operating  characteristics  of  practical  plasma  devices  such  as 
MUD  generators,  MPD  thrusters,  etc,  are  strongly  affected  by  several 
limiting  mechanisms  associated  with  current  transfer  between  electrodes 
and  the  plasma.  Particular  attention  should  be  drawn  to  the  excessive 
voltage  drop  and  power  dissipation  in  the  plasma-electrode  boundary 
layer,  and  the  electrode  erosion  associated  with  both  anodes  and 
cathodes . 

In  spite  of  intensive  research  over  many  decades,  the  plasma-  elec¬ 
trode  boundary  layer  is  still  poorly  understood  [2. 1,2.2].  The  diffi¬ 
culty  in  analyzing  this  region  stems  from  the  complexities  introduced  by 
the  interaction  of  electrical,  thermal,  fluid  dynamic  and  electrode 
surface  effects. 

This  work  is  directed  towards  providing  a  more  complete  under¬ 
standing  of  the  plasma-electrode  interaction  for  the  case  of  flowing, 
seeded,  thermal  plasmas  at  pressures  of  the  order  of  one  atmosphere,  and 
involves  both  theoretical  and  experimental  components. 

Consider  a  channel  flow  of  thermal  plasma  which  might  be  either  a 
noble  gas  or  combustion  products  gas  seeded  to  a  suitable  level  (~  1%) 
with  cesium  or  potassium.  We  may  define  the  electric  boundary  layer  as 
that  region  near  the  channel  walls,  composed  of  high  temperature  insula¬ 
tors  and  electrodes,  where  charged  particle  nonequilibrium  effects 
become  significant  and  even  dominant.  Normally,  the  walls  have  to  be 
operated  at  temperatures  less  than  that  of  the  core  flow  plasma  because 
of  material  limitations. 

On  physical  and  dimensional  grounds  the  region  between  the  core 
plasma  and  the  electrode  surface  can  be  divided  into  a  number  of  over¬ 
lapping  layers  of  decreasing  scale.  Proceeding  from  the  core  plasma  to 
the  electrode  these  layers  may  be  listed  as  follows  (see  Fig.  2.1): 


Figure  2.1.  Potential  distribution  across  plasma  diode. 

I  Equilibrium  Plasma 

II  Ionization  Non-Equilibrium  Layer 

III  Sheath.  (Not  to  scale) 

_2 

(i)  The  gasdynamic  boundary  layer  of  scale  6~10  m, 

depending  on  distance  along  the  channel,  and 

comprising  both  a  velocity  and  thermal  boundary  layer. 

(ii)  The  ionization  nonequilibrium  (or  ambipolar  diffusion) 

layer  of  scale  the  characteristic  ion 

recombination  length. 

(iii)  The  sheath  of  scale  equal  to  the  Debye  length  = 

(eokT/n0e2)1'  ,  typically  10-6m.  In  this  layer  charge 
neutrality  is  not  preserved  and  large  electric  fields 
develop. 

(iv)  The  Knudsen  layer  where  free  molecule  flow  takes 

over  from  a  continuum  (fluid)  description  and 
the  appropriate  scale  is  the  mean  free  path.  The 
mean  free  paths  for  electrons  and  ions  are  of  order 
£g  ~  ^i  ~  10”^®  and  are  thus  comparable  to  Xp. 
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Ic  should  be  emphasized  that  the  four  layers  distinguished  above  are 
overlapping  and  no  sharp  boundaries  can  be  distinguished.  Moreover,  the 
scale  lengths  are  nominal  ones  based  on  freestream  rather  than  local 
conditions.  In  general  the  scale  lengths  are  ordered  according  to  6  » 
tR  »  xD~  le  ~  i±. 

As  the  current  density  J  is  increased  from  zero,  the  potential  drop 
at  the  cathode  increases  and  more  ions  than  electrons  are  collected.  The 
reverse  happens  at  the  anode.  However,  these  conditions  only  hold  for 
very  small  currents  J  <  the  random  ion  current  (<  10  mA/cm^).  For 
higher  currents,  of  technical  interest,  the  cathode  must  emit  electrons. 
Of  the  known  emission  processes,  both  secondary  emission  due  to  ion 
bombardment  and  thermionic  (including  field-assisted,  Schottky  emission) 
are  relevant.  Depending  on  conditions,  particularly  the  cathode  mate¬ 
rial,  temperature  and  work  function,  diffuse  emission  persists  up  to 
some  current  density  when  there  is  a  transition  to  a  constricted  dis¬ 
charge  terminating  in  an  arc  cathode  spot. 

At  the  anode  it  might  be  expected  that  large  currents  up  to  J  ~  Jfi, 
the  random  electron  current  (~  10  A/cm  ) ,  could  be  collected  as  a  result 
of  the  reducing  sheath  potential  drop.  However  it  is  commonly  obse-.ec 
that  constrictions  in  the  form  of  anode  glow  spots  occur  at  lower 
currents  and  the  sheath  potential  drop  may  change  sign  as  indicated  in 
Fig.  2.1.  The  reason  for  this  behavior  is  unexplained,  but  it  is 
probably  related  to  the  creation  of  a  highly  conducting  discharge 
channel  which  penetrates  the  lower  conductivity  boundary  layer  plasma. 

From  the  foregoing  discussion  it  is  clear  that  the  physical  pro¬ 
cesses  Involved  are  extremely  complex  and  experimental  data  are  vital  to 
guide  the  theoretical  modeling. 


Research  Objectives 

1.  To  observe  and  record  the  current-voltage  characteristics  and  the 
transition  from  diffuse  to  constricted  discharge  modes  for  elec¬ 
trodes  immersed  in  a  flowing,  atmospheric  pressure  thermal  plasma. 
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2.  To  predict  the  current-voltage  characteristics  for  anodes  and 
cathodes  and  to  compare  these  predictions  with  the  experimental 
results. 

3.  To  use  the  Insight  provided  from  objectives  1  and  2  above  to  more 
fully  understand  the  constriction  process. 

Status  of  Research 

The  present  research  is  comprised  of  both  theoretical  and  experi¬ 
mental  components.  To  date,  the  theoretical  aspects  of  the  project  have 
been  stressed. 

Theoretical  Work 

The  approach  adopted  for  theoretical  modeling  is  to  start  off  with 
the  simplest  relevant  model  and  to  introduce  the  complexities  of  the 
complete  problem  one  at  a  time.  To  this  end,  the  plasma-sheath  problem 
for  a  plane  electrode  in  contact  with  a  uniform,  equilibrium,  weakly- 
ionized  plasma  has  been  formulated  in  some  generality.  By  assuming  that 
the  plasma  is  Isothermal  and  stationary,  the  complexities  associated 
with  the  gasdynamic  thermal  boundary  layer  and  Joule  heating  are,  for 
the  present,  eliminated  and  the  essential  plasma  physics  aspects  of  the 
problem  are  clearly  revealed. 

This  formulation  in  terms  of  the  conservation  equations  for  the 
electron  and  ion  concentrations  and  momenta  (including  inertia), 
together  with  Poisson's  equation,  gives  a  uniform  description  of  the 
ionization  nonequilibrium  layer  and  sheath,  in  both  the  colllonless  and 
collision  dominated  regimes,  without  having  to  resort  to  an  arbitrary 
division  and  matching  of  the  separate  regions. 

In  the  current  reporting  period  this  plasma-sheath  problem  has  been 
investigated  in  the  continuum  (collision-dominanted)  limit,  neglecting 
the  inertia  terms.  In  this  case,  the  governing  equations  are: 
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Here  x  is  normalized  on  the  ion  recombination  distance  ip  =  (2Da/f5n£)*^ 
where  0  is  the  ambipolar  diffusion  coefficient  and  6  is  the 

three-body  recombination  coefficient.  The  electron  and  ion  concentra- 
tions  are  normalized  on  n^  while  the  net  generation  rate  n  is  normal- 
ized  on  6  n*,.  The  potential  $  is  normalized  on  (kT/e)  while  the  elec- 

=  (kT/e£R).  The  fluxes  are  normalized  on 
rR  =  (2Dan<#/£R)  and  the  current  density  is  normalized  on  JgEen. 

(ye+^i) 

There  are  three  parameters  to  be  specified:  the  mobility  ratio 

y  =  (y1/ye)  «  1,  a  function  only  of  the  type  of  ion  and  background  gas; 

e  *  (xi./£r)  «  1.  the  ratio  of  the  Debye  length  to  the  ion  recombination 

length;  and  "J,  the  normalized  current  density. 
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where  a  is  the  electrical  conductivity  for  x  *». 


trie  field  is  normalized  on  E,? 


In  view  of  Eq.  (2.6),  four  independent  boundary  conditions  are 

required  for  solution  of  the  five  coupled,  first  order,  nonlinear  dif¬ 

ferential  equations  (2.1)-(2.5).  In  the  continuum  limit  two  of  these 
are  n  *  n  ■  0  at  r  •  0.  The  other  two,  E  and  either  r  or  F.,  are 

chosen  such  that  the  solution  goes  over,  at  large  x,  to  the  quasi¬ 

neutral  solution  derived  earlier  [2.3]. 


In  this  quasi-neutral  solution  we  put  n  ■  n,  ■ 

e  i 

Poisson's  equation  (2.5)  to  obtain  the  analytic  results: 
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^  -  in  tanh  x  -  T  in  sinh  x 
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It  is  noteworthy  that  If  has  a  simple  pole  and  $  has  a  logarithmic 

pole  at  x  -  0.  Moreover,  for  a  cathode,  for  3"  <  -  2y/(l  +  y),  T 

e 

becomes  positive  and  electron  emission  is  implied.  Furthermore,  for  an 
anode  ^  changes  sign  (corresponding  to  a  positive  anode  fall)  for 
J  >  (1  -  y)/(l  +  y),  and  F^  becomes  positive,  indicating  ion  emission 
for  T  >  2/(1  +  y). 

The  profiles  of  n*  and  n  are  shown  in  Fig.  2.2,  while  the  potential 
profiles  for  various  T  are  shown  in  Fig.  2.3. 

Integration  of  the  full  plasma-sheath  equations  (2.1)-(2.5)  has 
proven  rather  difficult.  Initially  a  "shooting"  technique  was  used,  in 
which  the  initial  values  of  E  (0)  and  F^(0)(or  1^(0))  were  adjusted  by 
iteration  to  obtain  the  solution  best  matching  the  quasi-neutral 


Figure  2.2.  Profiles  of  normalized  charged 
particle  concentrations  and  net  generation 
rate. 


i 

TIT/'r 


Figure  2.3.  Profiles  of  the  normalized  potential  as  a 
function  of  normalized  current  density:  7  <  0,  cathode 
7  >  0,  anode.  Broken  lines  indicate  the  asymptotes 
corresponding  to  the  uniform  electric  field  in  the 
distance  plasma. 


solution  (Eqs.  2.8-2.13)  at  large  x.  It  was  found  that  the  solution  was 
extremely  sensitive  to  the  choice  of  1T(0)  and  that  it  invariably  became 
unstable  for  x  >  2.5. 

Recently,  the  method  of  quasilinearization  with  orthonormalization 
has  been  utilized  whereby,  in  addition  to  the  boundary  conditions 
ne(0)  ■  0  ,  n^(0)  -  0  ,  the  boundary  conditions  are  imposed  that 
T^(x)  and  E  (x)  satisfy  Eqs.  (2.10)  and  (2.12)  at  some  large 
x,  e.g.  x  -  4. 

Typical  results  for  the  case  u  -  (1/300),  e  »  10_1  and  various 

values  of  7  are  shown  in  Figs.  2.4  through  2.7.  In  these  solutions  the 

quasi-neutral  ionization  nonequilibrium  region  merges  smoothly  into  the 

sheath  region  where  n^  *  ne.  One  can  arbitrarily  choose  to  take  the 

plasma-sheath  boundary  at  x^  where  In^  -  ne|/(ni  +  ng)  “6,  a  small 

parameter.  Referring  to  Fig.  2.4  for  7  -  0,  it  is  seen  that  a  reasonable 

choice  is  6  ■  0.1.  This  gives  a  sheath  thickness  of  x  -  0.82  or 

8 

xg  -  8.2Ad,  a  potential  fall  across  the  sheath  of  -  6.4  and  a  poten¬ 
tial  fall  through  the  ionization  nonequilibrium  layer  of  I  ■  1.1.  The 

K 

criterion  S  «  0.1  is  taken  as  corresponding  to  the  plasma-sheath 

boundary  in  all  the  cases  shown. 

In  Fig.  2.4  for  7-0,  it  is  seen  that  n.  and  n  lie  well  below 

i  e 

the  quasi-neutral  value  n"  -  tanh  x  out  as  far  as  x  «  3.  As  a  result 
the  region  of  net  generation,  n  moves  to  larger  x  than  in  the  quasi¬ 
neutral  solution.  Also  for  this  case  r  ■  r  everywhere,  tending  to  zero 

6  1 

at  large  x  and  to  -  1/2  at  x  ■  0.  The  net  generation  of  electron  ion 
pairs  generates  the  equal  fluxes  which  flow  to  the  electrode. 

Figures  2. 5-2. 6  show  solutions  for  anodes  (7  >  0)  as  the  current  is 

increased.  For  7  <  l,  the  potential  fall  in  the  anode  sheath  is  negative 

and  there  is  an  ion  sheath.  For  J  =  1,  illustrated  in  Fig.  2.5,  there  is 

no  sheath  and  n^  and  n^  closely  follow  the  quasi-neutral  approximation 

n  -  tanh  x  .  For  J  >  1  the  potential  fall  becomes  positive  and  there  is 

an  electron  sheath  as  is  clear  from  Fig.  2.6  for  7-2.  The  sheath 

thickness  and  potential  fall  for  this  case  are  x  -  0.47,  and  4  - 

8  8 

-1.87.  As  7  increases  the  electron  sheath  thickens  and  the  potential 
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2.6.  Solution  for  u  ■  1/300, 


drop  increases.  These  results  are  in  accord  with  the  results  from  the 
quasi-neutral  approximation.  All  values  of  the  potential  fall  through 
the  sheath  are  given  by  the  difference  between  the  calculated  potential 
and  the  potential  due  to  the  (assumed)  uniform  conductivity  of  the 
plasma . 

Figure  2.7  shows  the  solution  for  a  cathode  with  T  ■  -  2ji.  For  this 
case  the  electron  flux  to  the  cathode  is  reduced  to  zero  and  for  larger 
currents  1^(0)  becomes  positive,  indicating  electron  emission  as  pre¬ 
dicted  by  the  quasi-neutral  approximation.  For  T  ■  -  2 y ,  the  current  at 
the  cathode  is  carried  by  arriving  ions,  while  the  net  volume  generation 
of  electrons  provides  the  current  in  the  distant  plasma,  which  is 
carried  primarily  by  electrons.  Compared  with  the  case  T  >  0  (Fig.  2.4), 
the  cathode  with  T  ■  -  2ii  shows  an  increase  of  the  ion  sheath  thickness 
from  -  0.82  to  1.19  and  an  increase  of  sheath  potential  drop 

from  *  6.44  to  $  *  12.71.  With  larger  currents  the  cathode  sheath 

8  8 

thickness  and  potential  drop  can  be  expected  to  increase  significantly. 

Numerical  convergence  was  observed  to  be  quite  slow  for  anode  cases 
with  T  >  2  and  cathode  cases  with  T  <  -2y  .  This  has  led  to  a  re¬ 
examination  of  the  applicability  of  the  continuum  boundary  conditions, 
i.e.,  n^(0)  ■  n^(0)  “  0  .  It  is  worthy  to  note  that  the  convergence 

significantly  slowed  for  anode/cathode  cases  where  ion/electron  emission 
by  the  electrode  was  expected.  Considering  first  the  cathode  cases; 
electron  emission  may  be  expected  to  create  a  non-zero  value  for  n^  at 
the  electrode  surface.  Indeed,  by  constraining  only  n^(0)  »  0  and 

requiring  that  n  (x)  satisfy  Eq.  (2.8)  at  a  large  value  of  x  ,  the 
e 

calculations  predict  that  n  (0)  ■  0.091  for  T  ■  -2 p  .  The  resulting 

e 

solution  is  shown  in  Fig.  2.8.  This  implies  a  very  thin  electron  sheath 
immediately  adjacent  to  the  electrode  surface  followed  by  a  thicker  ion 
sheath.  Comparison  of  Figs.  2.7  and  2.8  shows  that  the  latter  solution 
follows  the  quasi-neutral  prediction  much  more  closely.  The 
sheath  thickness  and  potential  fall  are  found  to  be  "x  *  0.32 
and  ■  1.13  compared  with  the  much  larger  values  xs  ■  1.19  and 

$a  *  12.71  calculated  previously.  This  is  to  be  expected  because  now  the 
ion  sheath  near  the  electrode  has  been  partly  neutralized  by  electrons. 
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Figure  2.7.  Solution  for  n  ■  1/300,  e 
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It  is  clear  that  this  system  of  equation  is  extremely  sensitive  to  the 
electrode  surface  boundary  conditions. 

Application  of  the  same  boundary  conditions  for  the  T  “  0 
case  leads  to  the  solution  shown  in  Fig.  2.9.  Here  we  find  that 
x  “  0.31  and  J"  “1.13  compared  with  x  ■  0.ir<2  and  £  “  6.4  found 

S  S  S  8 

previously. 

As  the  current  increases  we  expect  the  sheath  to  disappear  near 

T  “  1.  Using  an  argument  similar  to  that  above,  we  choose  to  constrain 

ne(0)  “  0  and  require  n^(xj  to  satisfy  Eq.  2.8  at  a  large  value  of 

x  for  positive  currents  whereby  we  expect  ion  emission  from  the  anode. 

We  find  for  J  “  2  that  n,  (0)  -  0.071,  x  “  0.29,  and  ?  “  0.72.  This 

i  s  s 

may  be  compared  with  x  “0.47,  and  <b  “1.87  found  before.  This  solu- 

s  s 

tion  is  plotted  in  Fig.  2.10. 

The  above  results  show  the  critical  role  that  the  boundary  condi¬ 
tions  play  in  the  solution  of  this  problem.  The  derivation  of  physically 
realistic  boundary  conditions  when  there  is  electron  or  ion  emission 
from  the  electrode  represents  a  major  goal  of  this  project. 

Experimental  Work 


The  experimental  work  is  being  conducted  in  a  versatile  3  KW  thermal 
diffusion  flame  burner  capable  of  producing  a  combustion  plasma  to  3000K 
and  seeded  to  1%  by  mass  of  potassium  vapor  to  provide  a  conductivity  of 
up  to  10  S/m.  This  facility  was  built  in  connection  with  an  earlier  NSF 
grant  for  the  testing  of  plasma  diagnostic  techniques. 

The  burner  is  presently  constructed  with  a  water-cooled,  ceramic- 
lined  chimney.  The  chimney  inside  dimensions  measure  2.5  cm  by  8.75  cm. 
The  test  section  consists  of  a  small  volume  within  the  chimney  assembly. 
Optical  access  is  provided  by  two  2.5  cm  diameter  quartz  windows.  Gas 
temperature  measurements  throughout  the  test  volume  and  ceramic  surface 
temperature  measurements  have  been  made  for  a  variety  of  flow  condi¬ 
tions.  Typical  results  for  the  center-line  gas  temperature  and  ceramic 
surface  temperature  are  given  in  Figs.  2.11  and  2.12.  These  data  are  for 
a  relatively  cool  methane/air  flame  with  a  fuel/air  ratio  of  1.2. 
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Electrode  experiments  will  be  conducted  using  more  energetic  reac¬ 
tants  such  as  hydrogen/me thane /oxygen  mixtures  to  obtain  higher  temper¬ 
atures  and  electrical  conductivities.  The  center- line  gas 
temperature  for  a  portion  of  the  test  volume  (up  to  the  melting  point  of 
the  thermocouple)  and  the  ceramic  surface  temperature  for  such  a  flame 
are  plotted  in  Figs.  2.13  and  2.14.  These  data  were  taken  for  a  fuel 
mixture  with  a  hydrogen/ carbon  ratio  of  12  and  a  fuel/air  ratio  of  1.2. 

It  is  desirable  to  place  the  electrodes  approximately  4  -  5  cm  above 
the  burner  surface.  Unfortunately,  this  region  is  presently  marked  by  a 
steep  axial  gas  temperature  gradient.  For  this  reason,  an  improved  test 
section  is  being  constructed  with  smaller  optical  ports  and  increased 
thermal  insulation  of  the  chimney,  to  achieve  a  more  uniform  axial 
temperature  profile. 
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2.3  Interaction  of  Discharges  and  Pluid  Dynamics 
Introduction 

In  several  space  power  and  propulsion  systems  of  potential  interest 
to  the  Air  Force,  such  as  MFD  thrusters  and  MUD  generators,  fluid  mo¬ 
tions  induced  by  discharges  are  of  critical  importance  to  both  the 
understanding  and  performance  of  the  system.  In  any  MHD  device  the 
electromagnetic  field  may  affect  the  plasma  flow  field  in  several  ways. 
Among  these  effects  are  the  following: 

(1)  the  effect  of  Joule  heating  on  plasma  temperatures; 

(2)  the  Hartmann  effect; 

(3)  turbulence  suppression;  and 

(4)  the  generation  of  secondary  flows. 

The  first  three  of  these  effects  have  been  the  subjects  of  previous 
studies  performed  at  HTGL  [3. 1-3. 4;  these  are  summarized  in  Ref.  3.5]. 
The  fourth  effect,  secondary  flows,  is  the  area  in  which  our  current 
efforts  are  primarily  focused.  We  are  also  continuing  our  investiga¬ 
tions  in  the  area  of  turbulence  suppression  by  magnetic  fie'lds, 
including  the  effect  of  secondary  flows  on  turbulence  intensity  pro¬ 
files. 

The  physical  mechanism  for  secondary  flow  in  a  confined  plasma  flow 
with  a  transverse  magnetic  field  is  illustrated  in  Figure  3.1  For  sim¬ 
plicity  we  illustrate  a  distribution  of  the  Hall  current,  Jx(z),  which 
one  would  expect  from  a  first-order  analysis.  The  magnetic  field  B»B2 
is  assumed  constant.  The  resulting  distribution  of  the  J  x  B  Lorentz 
force,  acting  as  a  body  force  against  the  fluid  in  the  channel  cross¬ 
plane,  then  causes  a  secondary  flow  field  which  is  characterized  by 
large-scale  vortices.  It  is  believed  that  the  presence  of  this  sort  of 
transverse  flow-field  may  have  profound  implications  for  the  performance 
of  MHD  devices  [3.6-3.7J.  The  existence  of  secondary  flow  has  also  been 
inferred  from  data  obtained  at  HTGL  for  temperature  profiles,  voltage 
drops  and  heat  transfer  [ 3. 8 J • 


-  40  - 


Figure  3.1  Secondary  flow  mechanism 


Secondary  flow  and  che  Hartmann  effecC  are  closely  related,  in  that 
they  are  both  manifestations  of  the  appearance  of  the  Lorentz  force  in 
the  fluid  momentum  equation  for  a  collision-dominated  plasma, 

P  "’Vp-VT  +  JiI  (3.1) 

For  a  linear  MHD  device  B  is  assumed  to  point  in  the  z-dlrection;  thus 
the  Lorentz  force  has  components  in  the  x-  and  y-directions.  The  axial 
component,  Jyi,  gives  rise  to  the  Hartmann  effect,  while  the  transverse 
component,  -JXB,  causes  secondary  flow.  Whereas  the  Hartmann  effect  has 
been  of  Interest  since  Hartmann's  first  Investigation,  in  1937,  it  is 
only  recently  that  secondary  '  has  received  attention. 


The  presence  of  fluid  turuui  :e  alters  in  a  profound  way  the 
transfer  of  momentum,  energy  and  diffusion  in  the  fluid.  Thus,  any 


model  of  such  a  flow  must  incorporate  some,  necessarily  approximate, 
representation  of  turbulent  transport.  This  transport  is  inherently 
dependent  on  the  magnitude  of  the  turbulence  intensity. 

A  magnetic  field  acts  to  impart  Lorentzian  forces  on  the  motion  of 
the  conducting  fluid  in  general,  and  on  turbulent  motion  in  particular. 
These  forces  tend  to  suppress  or  damp  turbulent  components  of  velocity 
in  the  directions  perpendicular  to  the  magnetic  field.  Turbulence  damp¬ 
ing,  and  resultant  modifications  of  fluid  velocity  profiles,  are  well 
established  for  liquid  metal  flows  [3.9].  We  have  made  measurements  of 
the  turbulence  intensity  in  partially  ionized  plasmas  [3.4]  and  have 
found  a  similar  but  quantitatively  different  behavior.  These  results 
demonstrated  that  the  turbulence  intensity  can  be  reduced  by  approxi¬ 
mately  a  factor  of  two  for  the  moderate  values  of  magnetic  interaction 
achievable  in  a  continuous  combustion  MHD  device.  However,  it  is  not  yet 
clear  to  what  extent  this  turbulence  suppression  Influences  profiles  of 
mean  axial  velocity. 

The  problem  of  predicting  profiles  of  mean  axial  velocity  and  of 
turbulence  intensity  is  made  more  complicated  by  the  presence  of  dis¬ 
charge-induced  secondary  flows.  In  the  absence  of  secondary  flows  the 
effects  of  a  magnetic  field  on  turbulence  Intensity  and  mean  velocity 
profiles  are  distributed  symmetrically  over  the  channel  cross-plane.  The 
large  vortex  structures  characteristic  of  secondary  flow  may  cause  a 
convective  distortion  of  these  profiles.  For  example,  the  secondary  flow 
field  depicted  in  Figure  3.1  could  sweep  slow-moving,  turbulent  fluid 
into  the  top-wall  boundary  layer  while  pushing  fast-moving,  less  turbu¬ 
lent  core  fluid  downward  toward  the  bottom  wall. 

Furthermore,  the  effects  of  electromagnetic  forces  on  turbulence  in 
an  MHD  device  may  cause  the  turbulence  to  become  significantly  aniso¬ 
tropic,  because  Lorentz  forces  do  not  act  along  u.  Mc  field  lines. 
The  extent  and  distribution  of  such  anisotropies  are  im,  'rtant  factors 
in  turbulent  transport.  This  is  a  phenomenon  which  has  received  little 
experimental  attention. 
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Itmrch  Objectives 

The  over-riding  goal  of  our  work  la  this  area  is  to  obtain  experi¬ 
mental  information  which  will  allow  us  to  describe  the  three-dimensional 
flow  field  in  an  MUD  device.  Specifically,  the  major  objectives  which 
presently  guide  this  work  are  the  following: 

(1)  To  study  the  effect  of  electromagnetic  body  forces  on  the  secondary 
flow  field.  What  are  the  magnitudes  of  these  secondary  flows,  and 
how  is  their  vortex  structure  affected  by  applied  magnetic  fields 
and  discharges? 

(2)  To  study  the  effect  of  MHD  interaction  on  mean  axial  velocity 
profiles.  Are  they  affected  by  turbulence  suppression  caused  by  the 
magnetic  field,  and  to  what  extent  are  they  distorted  by  secondary 
flows? 

(3)  To  study  the  effect  of  MHD  interaction  on  turbulence.  What  is  the 
detailed  nature  of  magnetic  field  turbulence  suppression,  including 
the  development  of  anisotropic  turbulence?  To  what  extent  are 
turbulence  intensity  profiles  distorted  by  secondary  flows? 


Status  of  Kesearch 

Our  experiments  on  MHD  boundary  layers  and  secondary  flow  have  been 
performed  at  the  Stanford  M-2  MHD  facility.  The  flow  train,  shown  in 
Figure  3.2,  consists  of  a  2-MWth  combustor,  a  large  mixing  plenum  for 
the  combustion  products,  the  active  channel,  a  diffuser,  and  an  exhaust 
system  which  includes  a  scrubber.  The  active  channel  is  enclosed  in  a 
2.6-T  copper-coil  magnet.  The  top  and  bottom  walls  of  the  channel  can  be 
fitted  with  up  to  thirteen  pairs  of  electrodes.  These  are  connected 
externally  through  load  resistances,  and  can  be  configured  by  means  of  a 
patch  panel  so  that  various  voltages  are  applied  by  means  of  a  bank  of 
automobile  batteries.  This  allows  us  both  to  achieve  higher  levels  of 
magnetic  Interaction  and  to  gain  external  control  of  MHD  forces. 


As  our  experiments  on  secondary  flow  were  to  be  the  first  such 
measurements  in  an  MUD  channel,  it  was  desirable  to  design  experimental 
conditions  which  were  as  simple  and  well-controlled  as  feasible.  The 
most  straightforward  situation  is  the  one  suggested  In  Figure  3.1,  in 
which  the  net  Hall  current  is  non-zero.  Therefore  a  pair  of  electrodes 
at  the  upstream  end  of  the  channel  was  connected  externally  to  a  pair  at 
the  downstream  end,  as  shown  in  Figure  3.3.  The  magnetic  interaction  was 
then  maximized  by  driving  a  relatively  large  Hall  current  down  the 
length  of  the  channel  and  by  operating  at  a  low  flowrate. 
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Figure  3.2  M-2  facility  flow  train. 
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Figure  3.3  Electrical  configuration  for  the 
secondary  flow  experiments. 


Laser  Doppler  anemometry  was  used  to  measure  plasma  velocities  in 
both  the  x-  and  y-directions.  A  schematic  of  the  LDA  system  is  shown  in 
Figure  3.4;  a  more  detailed  description  is  found  in  Ref.  3.10.  Essen¬ 
tially,  a  conventional  single-component,  dual-beam  system  was  adapted  to 
the  particularly  difficult  circumstances  of  a  combustion  MHD  flow. 

Preliminary  measurements  of  secondary  flow,  made  in  1983,  are  dis¬ 
cussed  in  Ref.  3.10.  These  were  made  under  a  variety  of  conditions,  but 
were  limited  to  one  quadrant  of  the  channel  cross-plane.  During  the  past 
year  we  have  made  more  detailed  and  comprehensive  measurements , 
including  measurements  of  the  y-directed  velocity  component  over  the 
entire  channel  cross-plane.  With  all  other  conditions  constant,  the 
magnetic  Induction  was  varied  to  produce  three  cases:  (1)  B  -  2.4  T, 
pointing  in  the  positive-z  direction;  (2)  B  “  2.4  T,  pointing  in  the 
negative-z  direction;  and  (3)  B  *  0.  The  corresponding  effect  on  the  J  x 
B  body  force  experienced  by  the  plasma  in  the  cross-plane  is  that  (1)  it 
pointed  in  the  negative-y  direction;  (2)  it  pointed  in  the  positive-y 
direction;  and  (3)  it  was  zero. 

Figure  3.5  shows  the  results  for  the  positive-B  case.  The  picture 
here  appears  immediately  to  confirm  that  Lorentz  forces  strongly  affect 
the  transverse  flow  field.  The  most  significant  qualitative  feature  of 
the  flow  field  is  the  degree  to  which  it  conforms  to  the  first-order 
model  of  Figure  3.1:  the  top  part  of  the  cross-plane  is  characterized  by 
two  large,  counter-rotating  vortices,  which  in  the  core  follow  the 
negative-y  direction  of  the  Lorentz  force.  The  bottom  of  these  vortices 
(where  the  secondary  flow  is  directed  in  the  z-direction,  and  hence  the 
y-component  is  zero)  is  located  at  Y/H  2  0.3.  The  peak  measured 
velocity  had  a  value  12%  of  the  bulk  axial  velocity.  This  velocity 
evidently  peaks  quite  close  to  the  sidewall,  so  that  the  actual  maximum 
was  undoubtedly  larger  than  this;  in  the  negative-B  case,  in  which 
measurements  were  made  closer  to  the  sidewall,  the  maximum  measured 
value  was  14%.  The  region  near  the  bottom  wall  does  not  appear  to  be 
characterized  by  a  vortex  rotation,  but  rather  by  a  lifting  of  the  flow 
away  from  the  wall.  This  behavior  is  not  predicted  by  the  first-order 
model,  and  warrants  further  study. 


SCANNING 

MIRROR 


BEAMSPLITTER 


Figure  3.4  Laser  doppler  anemometry  system. 

Measurements  made  on  completely  separate  occasions  (but  under  the 
same  operating  conditions)  have  shown  that  the  data  were  quite  reproduc¬ 
ible.  Furthermore,  measurements  made  in  the  negative-B  case  showed  the 
same  results,  but  rotated  by  180°  about  the  channel  axis;  this  is  pre¬ 
cisely  what  we  expect,  as  reversing  the  z-direction  of  the  magnetic 
field  reverses  the  sign  of  the  y-component  of  the  Lorentz  force.  Final¬ 
ly,  the  zero-B  results  showed  a  much  weaker  secondary  flow  field,  with 
peak  velocities  5%  of  the  bulk  velocity,  and  without  the  two-cell  struc¬ 
ture  of  the  B-field  cases.  Taking  the  results  for  the  three  cases 
together,  the  evidence  appears  quite  strong  that  Lorentz  forces  can 
exert  a  definite  and  strong  influence  on  the  secondary  flow  field  in  an 
MHD  channel. 

We  would  expect  the  presence  of  such  a  strong  secondary  flow  to  have 
major  effects  on  the  plasma  thermal  and  momentum  behavior.  For  example. 
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the  vortices  seen  la  Figure  3.5  would  sweep  relatively  cold  sidewall 
fluid  into  the  top  wall  boundary  layer  while  pushing  hot  core  fluid 
downward  toward  the  bottom  wall.  Preliminary  inspection  of  the  data  for 
the  three  cases  appears  to  demonstrate  that  this  convective  effect  of 
secondary  flow  does  indeed  cause  a  redistribution  of  the  transverse 
profiles  of  axial  mean  velocity  and  turbulence  intensity. 
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Figure  3.5  Measurements  of  the  y-directed  velocity  component 
in  the  positive-B  case. 
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Figure  3.6  shows  measured  mean  axial  velocities  as  a  function  of  y 
along  the  z-centerline.  The  skewing  of  the  profiles  in  the  B-field  cases 
is  quite  pronounced.  Extending  these  measurements  to  the  y-z  cross¬ 
plane,  an  array  of  data  for  the  positive-B  case  yields  the  isovel  graph 
shown  in  Figure  3.7.  We  see  that  peak  axial  velocities  are  pushed  down¬ 
ward  toward  the  bottom  wall,  while  slow-moving  fluid  is  swept  over  the 
center  of  the  top-wall  boundary  layer  and  then  pushed  downward. 

If  the  fluid  being  swept  from  the  sidewall  into  the  top  wall  bounda¬ 
ry  layer  (positive-B  case)  is  cold  and  electrically  resistive,  it  is 
also  relatively  dense,  slow  and  turbulent.  Again,  our  preliminary 
analysis  of  measurements ,  with  conditions  corresponding  to  the  three 
secondary  flow  cases,  supports  our  expectations  in  this  regard.  Our  data 
indicate  that  the  region  of  minimum  freestream  turbulence  is  shifted 
downward  from  the  center  of  the  cross-plane  toward  the  bottom  wall;  its 
location  corresponds  quite  closely  to  the  shifted  location  of  the  axial 
velocity  peak. 

In  conclusion,  the  past  year  has  been  one  of  significant  achievement 
for  our  program  in  MHD-induced  secondary  flow.  Previous  to  this  program 
MHD  secondary  flow  was  merely  a  prediction.  Our  experimental  results  are 
showing  that  this  is  Indeed  a  significant  phenomenon.  The  data  from 
these  experiments  are  still  being  analyzed,  but  it  is  already  evident 
that  several  of  the  effects  of  secondary  flow  (on  profiles  of  axial 
velocity,  turbulence  intensity  and  electrical  conductivity)  have  been 
demonstrated  to  be  quite  pronounced. 
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